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Abstract. We employ the Q representation to study the non-classical correlations that are present from 
below to above-threshold in the degenerate optical parametric oscillator. Our study shows that such corre- 
lations are present just above threshold, in the regime in which stripe patterns are formed, but that they 
also persist further above threshold in the presence of spatially disordered structures. 

PACS. 42.50.Lc Quantum fluctuations, quantum noise, and quantum jumps - 42.50.Dv Nonclassical field 
states; squeezed, antibunched, and sub-Poissonian states; operational definitions of the phase of the field; 
phase measurements - 42.65.Sf Dynamics of nonlinear optical systems; optical instabilities, optical chaos 
and complexity, and optical spatio-temporal dynamics 



1 Introduction 

Quantum correlations in transverse spatial patterns of op- 
tical systems have been much studied in recent years. This 
is because of their fundamental relevance as a macroscopic 
quantum phenomena in spatially extended systems, and 
also because of possible applications for quantum image 
processing In this context a prototype system 

studied is the Optical Parametric Oscillator (OPO). Pat- 
tern formation in OPO was predicted long ago [pi, and 
there are recent experimental observations of spatial mode 
interaction and pattern formation JgJ . 

Quantum correlations in the OPO are often under- 
stood in terms of the bàsic process of emission of twin sig- 
nal photons in parametric down conversion of a pump pho- 
ton Below, but close to the threshold for paramet- 

ric oscillation, any critical transverse mode k c is weakly 
damped. The ring with radius |k c | seen in the far field is 
a noisy precursor Q of the wave-number to be selected 
above threshold. Spatial correlations have been studied in 
this regime showing interesting non-classical behavior [11|. 
This is referred in the literature as a quantum image [1C|. 
In particular, the intensity difference between any two op- 
posite modes with radius |k c | in the far field is subpoisso- 
nian. This is onc manifcstation of quantum cntanglement 
of the modes |ï^| . From a technical point of view, corre- 
lations in the regime below threshold can be calculated 
analytically within linear approximations and using the 
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Wigner representation [Q. Above, but relatively close to 
threshold, a stripe pattern appears (Sect. [?]) due to the in- 
terference of two selected signal modes of opposite critical 
transverse wavenumber p[ . Again, this result can be un- 
derstood as a macroscopic consequence of the microscopic 
process of twin photon emission. Close to threshold this 
phcnomenon is studied with models which consider only 
a small number of modes; only the pump mode and the 
critical ones are considered. The calculation is based on 
a linearization approximation and the use of the Wigner 
representation. At threshold or extremely close to thresh- 
old, critical fluctuations would invalidate a linearized cal- 
culation, but a few mode approximation, such as those 
considered in p~5| | , could be appropriate. 

Quitc a diffcrent situation exists when multimode non- 
linear dynamics have to be considered. Some results have 
been obtained in the context of multimode interaction in 
linear approximations in other systems |Ï^,|Ï^,|Ï8|. How- 
ever, a full multimode nonlinear description is needed well 
above threshold where the emission of the critical modes 
necessarily stimulates spatial harmonies in both the pump 
and the signal, and undamped modes associated to spatial 
broken symmetries play an important role. This is also the 
case in the convective regime caused by walk-off, where 
a noise sustained pattern appears with a broad (multi- 
mode) spectrum and amplified nonlinear fluctuations [ fÏ9[ 
pfï| . In these situations two relevant qüestions have to be 
addressed. First, from the conceptual point of view one 
might ask if the quantum correlations are degraded by the 
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subsidiary nonlincar processes of multimode competition. 
In particular, can the twin beams that appear close to 
thrcshold bc incoherently depleted by other processes? If, 
however, quantum effects persist, then the simple explana- 
tion in terms of emission of pairs of twin photons needs to 
be revised. Secondly, and from a technical point of view, 
new calculation techniques or approximations that go be- 
yond linearized approximations need to be developed. 

We have previously examined these qüestions using 
a time-dependent parametric approximation |2^,^T|. We 
concluded that quantum correlations were strongly sup- 
pressed in the noise-dominated convective regime ]2(J . How- 
ever, in the regime of absolute instability we found that 
non-classical correlations between the ±fc c signal beams 
persists in spite of nonlinear interactions with the pump 
and with higher order harmonies of the signal Q] • This is 
still true in cases in which the mean intensity of the two 
critical signal beams is different because of walk-off [gïj . 
Two limitations of our time dependent parametric approx- 
imation, which is well suited in the convective regime, are 
that fluctuations in the pump are neglected and that its 
validity is restricted to vàlues relatively close to threshold. 

As an alternative approach we propose, in this paper, 
the use of the Q-representation and its associated nonlin- 
ear Langevin equations for c-number complex fields. With 
this method we can study quantum spatial correlations in 
the DOPO without any linearization, few mode approxi- 
mation or time-dependent parametric approximation, for 
vàlues of the pump less than twice its threshold value. In 
particular we can reach high pump vàlues for which the 
stationary classical solution for the signal field is not a 
stripe pattern, but rather a homogeneous state with one 
of two preferred phases. Spatial coexistence of domains 
of these two homogeneous solutions separated by domain 
walls gives rise to spatially disordered patterns (seemingly 
chaotic). These have a broad spectrum in the far field for 
both the pump and the signal field. This structure can- 
not be described, not even as a first approximation, as 
the interference of two twin beams of opposite transverse 
critical wavenumber. This is because of other cascading 
processes coupling many different modes and freqüències. 
Nevertheless, we find that any ±fc pair of modes in the 
far field that lie within the active range of the broad spec- 
trum show non-classical correlations. This seems to indi- 
cate that the bàsic process of parametric down conversion 
is behind correlations found even in these very compli- 
cated spatial structures. 

The paper is organized as follows. The quantum for- 
mulation of the problem is presented in Sect. |2|, w here 
phase space methods are briefiy compared (Sect. |2.lj) and 
Langevin eq uatio ns for the Q representation are intro- 
duced (Sect. In Sect. | we briefiy review relevant 

classical results concerning the instabilities to pattern for- 
mation and to homogeneous solutions well above thresh- 
old. Our general aim is to study the quantum properties of 
the correlations in regimes that have not been previously 
studied. These include the critical point in presence of a 
multimode interaction (Sect. [|), and the above threshold 
region in presence of complex patterns (Sect. FJ). We find 



signatures of entanglement of the beams in the presence 
of stri pes influenced by the presence of higher harmonies 
(Sect. ^j]). We show that such entanglement is not only 
associated with a stripe pattern, but it i s present also in 
spatially disordered structures (Sect. 5.2). Finally Sect. || 
is devoted to concluding remar ks. 



2 Quantum formulation of DOPO dynamics 

To describe the intracavity dynamics in a DOPO we in- 
troduce the boson spatial modes ^4o( x , t) and ^4i(x, í), re- 
spectively at the pump frequeney 2u>, and signal frequeney 
lo, and satisfying Standard equal-time commutation rela- 
tions íl 



iiCx^iVx',*) = <5y<5(x-x'), i,j = 0,1 (1) 



Here x denotes the transverse coordinate(s). A Hamil- 
tonian operator describes the interaction between these 
modes in the non-linear medium. The intracavity fields 
constitute an open device (2^,|2í|, modeled within a sta- 
tistical approach in t he Schròdinger picture by a Master 
equation. In Sect. 2.1 we review the Master equation of a 



DOPO. We then report on the possible phase-space de- 
scriptions, introducing the Q-rep res entation and the asso- 
ciated Langevin equations (Sect.^1). 



2.1 Master equation and phase space descriptions 

The intracavity dynamics of our open system is described 
by a Master equation for the reduced density operator p 



(2) 



We consider a plane one-sided cavity, hence the Liouvillian 
accounting for dissipation through the partially reflecting 
mirror is given by 

^ — n i " 



The Hamiltonian operator, expressed as a function of fields 
operators Aq(x, t) and Ai(x, i), is: 



H = H + H int + H ext 



where [24 



H a = h í d 2 xJ2 [ldt(x)(A - a,V 2 )i,(x) 

describes free propagation of fields in the cavity, 
H ext = iH í d 2 *E fij(x) - i (x) 



(3) 



(4) 



(5) 
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is due to the interaction with the external pump E, which 
we choosc to be real, and 



Hint = ih 



dh 



Í0(x)ií 2 (x)-4 (x) i2 (x) j (g) 



is the interaction term between first and second harmònic. 

Dcnsity operators in state space can be mapped to 
quasi-probability distribution densities on phase space. 
These can be used to calculate ensemble averages of opera- 
tors in defined orderings p5| , p6|^ Using this "quantum-to- 
classical" correspondence, Eq7(|2[) can be converted into an 
equation of motion for a quasi-probability distribution in 
the phase-space of radiation fields «i(x), associated with 
the operators A;(x). 

The presence of non-linearities leads to a functional 
differential equation for the quasi-probability that is not 
of the Fokker-Planck type, leading to difficultics in ob- 
taining solutions J2?j. The offending term for the system 
of interest is Hi n t, which gives a functional term of the 
form 



[io(x)i{ (x)-h.c.,p\ 



sao- 



(7) 



Here s denotes the operator ordering selected. We see that 
third order derivatives appear in the temporal evolution of 
the Wigner representation (s = 0). The approximate equa- 
tions obtained simply dropping these third order terms 
constitute the basis of stochastic electrodynamics ]2§| ] . This 
approach works well in linear regimes, in which the Wigner 
distribution satisfies a genuine Fokker-Planck equation. In 
particular, in the DOPO below the threshold of signal gen- 
eration the intensity of the signal is of the order of the 
quantum noise, while the pump has a macroscopic mean 
value, so that its fluctuations can be neglected. With the 
assumption of a classical undepleted pump, the Hamilto- 
nian that describes the quantum dynamics of the signal 
is quadratic and the Langevin equation |27|] - equivalent 
to the Fokker Planck equation for Wigner representation 
- can be analytically solved |Ï3[ |. Recent investigations 
have shown the limits of this stochastic electrodynamics 
in reproducing quantum higher-order moments |29j| . The 
same type of approximation is generally possible above 
threshold, linearizing around a pattern solution fïl| . Due 
to the in-homogeneity of the reference state, moreover, 
only semi-analytical or numerical simulations can be pro- 
vided. Above threshold great care has to be taken if the 
system is translational invariant (flat mirrors and homo- 
geneous transversal pump profile). The pattern solution 
breaks the translational symmetry and therefore there is 
a Goldstone mode which is neutrally stable [^6| . Noise ex- 
cites this mode, giving diffusion of the phase which fixes 
the position of the pattern Jï(| . Moments involving such 
big fluctuations cannot be correctly described within a lin- 
earized treatment in the fields amplitudes. In particular, 
such an approach leads to unphysically divergent quadra- 
ture correlations, although correct results can be obtained 
for the intensity correlations [fl8| . 



Another well-known quasi-probability is the P+ repre- 
sentation 1 30 , consisting in the extension of the normal 
ordered P representation over a doubled phase space [^TJ . 
The P representation (s = 1 in Eq. (0)) suffers of nega- 
tive diffusion in problems of interest, but the P+ generally 
gives good results, with the advantage of the possibility to 
obtain immediately also the moments outside the cavity. 
However, in some systems this doubling phase-space tech- 
niquc has shown divergent trajectòries, as reviewed in Rcf. 
[ fj2[ . In particular, there are regimes in extended systems 
-like the convective regime [EOl — in which the presence 
of large fluctuations around the unstable reference state 
would result in diverging trajectòries in the P+ and alter- 
native methods are needed. 

In this paper we employ the Q representation corre- 
sponding to anti-normal ordering of field operators. The 
most important property of this representation is that it 
satisfies the requirements for a true probability distribu- 
tion. In fact the Q-representation may be defined as the 
diagonal matrix elements of the density operator in the 
space of coherent states 



Q{a ,ai) = — < a , ai|/5|a , ai > 

7T 



(8) 



and so is both positive and bounded |2(| . Due to the over 
completeness of the coherent states ensemble, the defini- 
tion (||) uniquely determines the density operator p. Phys- 
ically this representation, resulting from a Gaussian con- 
volution of the Wigner representation, corresponds to si- 
multaneous measurements of orthogonal quadratures, as 
limited by the Heisenberg principle, in a eight-port homo- 
dyne detector j33j. From Eq. (^|) (with s = — 1) we observe 
that the Q-representation suffers of negative diffusion. Un- 
like the Wigner function, however, the Q function is always 
positive and well-behaved. The possibility to obtain a pos- 
itive solution in presence of a negative diffusion [|4) lies 
in the presence of a restricted ensemble of initial condi- 
tions, that cannot be arbitrarily narrow. In other words, 
not all mathematical forms for the Q function correspond 
to physical states. The evolution of a physical state - cor- 
responding to an hermitian density operator - in the Q 
representation will always be positive [^6||3^]. We should 
note that a Q representation with a doubled phase-space 
has been proposed in order to deal with negative diffusion 
|p6| . This has been shown togive good results in some 
non- linear quantum systems [p7| . 

In the next section we investigate the possibility to 
use the Q-representation for devices consisting on a cav- 
ity filled with a x 2 medium, as in the OPO and Second 
Harmònic Generation (SHG). 



2.2 Langevin equations in Q-representation 

The evolution equation of the functional Q for our model 
(described above) is: 



dQ(ao,ai) 
dt 



d 2 x< - 



S 5 

-V a + - — —Vi+c.c. 
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S 2 S 2 
7o ía (x)(5a*(x') + 71 Sa^Sa*^')' 
5 2 

g&o- — r-77 — r-TT + c - c 



condition (|Ï4|), should these occur. £o is a white Gaussian 
noise with non-vanishing moment: 



where the drift terms are 



Q(ao, «i), 



Vq = -7o[(l + iA))-ïaoV 2 ]a (x,i)- ^a 2 {*,t) + E 

Ví = -71 [(1 + iAi) - ia 1 V 2 ]ai(x,í) +ga (x, í)aí(x,í). 

If the difïusion term is positive then our evolution equa- 
tion is a bona-fi.de Fokker-Planck equation. In the other 
case this equation doesn't describe an ordinary difïusion 
process. For equation (^) the difïusion term is positive ií 



0) 



(£ (x,íX *(x^')>=5(x-x')(í-í') 



(15) 



\a (x,t)\ < 



271 
9 



The modulus oí the stationary field at threshold takes the 
value \A t í ^ Lr \ = 71/.9 (see Sect. [|). This means that the 
condition (nG) corresponds to pump trajectòries taking 
vàlues that are less than twice the threshold value. Stay- 
ing in a region far from the limit (|ï(ï| ) - we are considering 
Aq* < l.Svlg - an extremely large fluctuation in a trajec- 
tory would be necessary in order to lose the positiveness oí 
the difïusion. Clearly these trajectòries have a negligible 
probability to appear, and never appeared in our simula- 
tions. For these reasons, the approximation we propose is 
to study Langevin equations related to the Fokker-Planck 
equation given by (^) and (|ïo|), neglecting any trajectò- 
ries that would make negative the difïusion term. Clearly 
the condition (FÍO|) does not depend on the írequency at 
which the system is pumped. For this reason the method 
is suitable, and has been already successfully used, to de- 
scribe non-linear fluctuations in stripe patterns in SHG in 
a regime of pump vàlues limited by Eq. (|ï^) pq ], 

From Eqs. (|s)-|lO|) , with the scaling (D is the transversal 
dimcnsionality oí the system) 

7o=7i=7 ; a o = «i/2 = a, (11) 

/ / x 

t = 7Í , x = —=, 
y/a 

ai 9_ à p/ _ 9_ rp 1 _ 9 

we obtain the equations: 

S t a (x, t)=- [(1 + iA ) - ÍV 2 } a (x, t) + E - 

5"i(x,í) + J--ío(x,<) 
2 V a 7 

<9 t ai(x,i) = - [(1 + Í4) - 2iV 2 ] ai(x,í) 

2 3, 



(12) 



ao(x,íK(x,í) + V--a(x,í). (13) 
V a 7 

The condition (|ï^) in the new variables is 

K(x,í)|<2. (14) 

We solve these Langevin equations by numerical simula- 
tion, neglecting any trajectòries that do not satisfy the 



The signal noise Çi results to be phase sensitive, due to 
the presence oí diagonal terms in the difïusion matrix of 
Eq. (|^). Moreover this noise is multiplicative, depending 
on the value oí the pump field. However due to the form 
of H (quadratic in A\ and linear in Aq), these equations 
have the same formal cxpression in the Ito or Stratonovich 
interpretations J27J. 

The phase sensitive multiplicative noise £i(x, t) can be 
written as 



( 10 ) £i(M) = 



-a 0I (x,t) 



+-^/2 + a R(x,t) 

2^2 + a 0R (x,t) 2 




\ao(x,t)\ 2 



oíor(x,í) 



ip{x,t) 



<j)(x,t) 
(16) 



with ao = aoR + iaai and 0,-0 uncorrelated real white 
noises in space and time, with variances one. Actually, the 
difïusion matrix of a Fokker-Plank equation fix only three 
of the four degrees of freedom in the choice of the real 
and imaginary parts of the noise term £1 ; this depends on 
the multiplicity of Langevin processes associated with the 
same Fokker-Planck equation ^7j. Hence Eq. (^6|) corre- 
sponds to a particular choice among several possible rep- 
resentations p9[ . 

In the next sections we present the numerical results 
obtained by simulating Eq. (jï^) and Eq. (|Ï3|) with the 
same integration method of RcL· jÏ6,38|. In particular we 
consider one transversal dimcnsion [D = 1) and parame- 
ters: 



Ao - 0, A 1 



-0.18, 



10" 



(17) 



with a system size of four critical wavelengths L = 4A C , 
with X c = y- (see Sect. |). 

Eqs. (12-Ï3|) are suitable for the study of quantum 
fluctuations in different regimes - discussed in Sect. ^ - 
from the linear regime below threshold to the multimode 
regimes quite above threshold. In Sect. |] we present cal- 
culations of quantum correlations below threshold and at 
the critical point, comparing with analytical results in the 
linear approximation. In Sect. |^ we consider the above 
threshold regime, in the situation of regular stripe pat- 
tern formation and also when disordered structures are 
formed (see Fig. pi). 



3 Pattern formation in OPO 

We first review the patterns that arise in the OPO in the 
classical limit. The phase matched DOPO, where both 
pump A and signal A\ fields are resonated, is described 
classically by equations identical to Eqs. ( |Ï2|) and (^) but 
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neglecting the noise terms. The classical equations have a 
trivial homogeneous solution 



Af = 0, 



ASt _ 

A Q — 



E 



1 + iA 



(18) 



A linear stability analysis around this solution gives the 
following dispersion relation for the growth of the signal 
ficld perturbations with wave vector k (the pump field Aq 
is always stable) [|j 



Ai(k) 



-1 ± 



K l \ 2 



{Ai + 2k 2 



(19) 



For negative signal detunings, the zero homogeneous so- 
lution becomes unstable at E = E c = ^\ + Aq. The 
perturbations with maximum growth rate are those with 
wave number \k c \ = ^/—Ai/2, and a pattern with this 
wave number is formed at threshold For positive sig- 
nal detunings the zero homogeneous solution is stable for 
E < E c = + Al)(l + Af). In this case the instabil- 
ity takes place at zero wave number leading to a non zero 
homogeneous solution [EÜI 



Af = ±. 



E\A{ 



st 12 



ASt 



t (l + i^o)(l+*^i) 
Ejí + iAj) 
(l+^o)(l + ^i) + |Af| 2 



\Af | 2 



(20) 
(21) 



There are two equivalent solutions for the signal held with 
a 7r phase difference. 

In this paper we will consider the case of zero pump 
detuning (A = 0) and negative signal detuning (Z\i < 0), 
in which stripe pattern arises at threshold (E c = 1). In- 
creasing further the pump the nonzero homogeneous so- 
lutions (^o|) become stable. In fact, we observe they are 
stable for pump vàlues around E = 1.1. This is lower 
than the valuc at which the stripe pattern becomes lin- 
early unstable jy],[Ï2|. Numerical studies in this regime 
find multistability between the stripe pattern, the two 
nonzero homogeneous solutions and several irregular spa- 
tially modulated solutions. The latter are formed by fronts 
with oscillatory tails connecting the two equivalent ho- 
mogeneous solutions. In systems with two spatial dimen- 
sions, there are also coexisting labyrinthine patterns. In 
this case (D = 2) both the irregular spatially modulated 
solutions and the labyrinthine patterns cease to exist at 
the modulational instability of a flat front connecting the 
two homogeneous solutions |^,^3). Here we consider sys- 
tems with only one spatial dimension for which such an 
instability does not exist. This means that the regime of 
multistability in parameter space is much larger. The ir- 
regular spatially modulated solutions found in this system 
are an example of frozen chaos as described in Ref . |Q . 
In that case the interaction of two distant fronts can be 
described by a potential with several wells which become 
progressively deeper as the distances between the fronts 
decreases. The OPO cannot be described in terms of such 
potentials, but numerical studies reveal equilibrium dis- 
tances whenever the màxima (or the minima) of the local 
oscillations of the front overlap with each other 



4 Quantum correlations below and at 
threshold 



The spatiotemporal dynamics of the signal field is shown 
in Fig. |ï] for two relevant vàlues of the pump, below but 
near to threshold (quantum images regime), E = 0.999 
and at threshold (E = 1). The far field (FF) shows strong 
f luctuat ions dominated by the critical wave- vector: in Sects. 
(11-12) we discuss the quadratures and intensity quan- 
tum correlations of these modes. 



4.1 Quadrature correlations 

The direction in which quadrature squeezing appears is 
determined by the eigcnfunction V±(k,—k) of the linear 
problem dtV±(k, — k) = X±(k)V±(k,—k), as reported in 
Ref. M: 



V±(k,-k) 



e iK ± {k) = T 



--e i0± 5A 1 (k)±5A* 1 (-k) 

iAi + 2ik 2 t VPgF ~ (^i + 2fc 2 

^0 



(22) 



The solution V+(k,—k) gives the direction of amplifica- 
tion of fluctuations, while fluctuations are damped for 
V-(k, —k), giving rise to quadrature squeezing. In partic- 
ular, for the critical wave-vector k c and for our choice of 
parameter (real Af^) we obtain V±(k c ,—k c ) = 5Ai(k c ) ± 
5A\(—k c ). Therefore, the largest squeezing at threshold 
will be in the difference of real parts and the sum of imag- 
inary parts of the ficld for wave-numbers k c and — k c . 
We define the real quadrature operator: 



X(k) =i 1 (fc) + i t 1 (fc) 

and the quadrature superpositions 

X-(k) = X(k)-X(-k) 
X+(k) =X(k)+X(-k), 



(23) 

(24) 
(25) 



corresponding, respectively, to damped and undamped quan- 
tities at threshold. 

Below threshold, within a linearization approximation 
[ p"3| , the normal-ordered variances normalized to the shot 
noise (N'x) 



arc: 



<: (X-(fc c )) 2 :> 
<: (X + {k c )f :> 



A/j 



x 



— E 
1 + E 

E 
1-E' 



(26) 
(27) 



These quantities coincide with the variances since the mean 
vàlues are zero: < X±(k) >= 0. The normal ordering al- 
lows us to immediately identify non-classical features asso- 
ciated with squeezing such as negative variances. Eq. ( |2^ ) 
shows an increasing degree of squeezing, approaching the 
value —0.5 at threshold. In Fig. || theoretical predictions 
and numerical results are shown to be in good agreement, 
confirming the validity of Eqs. ( |Ï2] - |T^ ) below threshold. 
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On the other hand Eq. ( p7| ) is always positive indicating 
that thc thc fluctuations in the direction of instability are 
essentially classical and larger than those found for a co- 
herent state. In Fig. || we show the agreement between 
theoretical predictions and numerical results for the un- 
damped quadrature, even as close as 1 %oto thrcshold. 
The limits of the linear treatment, discussed above, are 
now evident in the divergence of Eq. (|2^) for E — » 1 . In 
contrast, numerical simulation of the nonlinear Eqs. ([Ï2] - 
|Ï3| ) gives the expected saturation at the critical point, at 
a value which depends on the noise level. 



4.2 Intensity correlations 

We can find non-classical features in the intensities of the 
twin beams by evaluating the normal-ordered variance in 
the difference of the two intensities: 



V(k) 



[SN^k)- SN^-k)} 2 



(28) 



normalized to the corresponding shot noise value A/jv(fc). 
This value is proportional to the sum of the intensities of 
the two beams with wavevectors ±fe. Negative vàlues of 
V indicate sub-Poissonian statistics for thc intensity dif- 
ference of the two signal beams at ±k pïj. In a linear 
analytical treatment below threshold V(k) = —0.5, inde- 
pcn dcntly of the pump intensity and of the wave-vector 
|l3|,pT|. In other words the normalized intensity correla- 
tions, Eq. (p8|), do not show a non-classical behavior which 
is stronger for the critical wave vector or at the critical 
point. This is in contrast with the behavior of the quadra- 
tures correlations Eqs. ( p6||2"7| ). Nevertheless, the critical 
conditions are of significant interest because of presence 
of highcr intensities. 

The numerical expression of V(k) for different spatial 
modes (0 < k < 5fc c ) is compared, in Fig. ||, with the 
analytical value —0.5 for pump E — 0.99, showing good 
agreement. Small deviations for large wavevectors k can 
appear numerically due to the smallncss of the shot noise 
to which the variance is normalized; A/jv is proportional 
to thc mcan intensity of the field, shown also in Fig. || 

In Fig. H we plot the variance V(k c ) at the critical 
wave-vector as a function of thc pump E: Wc obtain good 
agreement with analytical predictions below threshold. 



5 Quantum correlations above threshold 

The non-linear equations in the Q-representation are used 
here to study the regime of pattern formation above thresh- 
old. These equations improve the time- dep endent para- 
metric approximation introduced in Ref. |2Cfl , which is well 
suited to study the convective regime (in presence of walk- 
off) or the instability point, because the pump fluctuations 
are disregarded with respect to the signal ones. Eqs.(|Ï2[) 
and ( |Ï3| ) are vàlid in a wide region above threshold, and 
give a complete description of the pump fluctuations. 



Slightly above the thrcshold (E = 1.02) the variance of 
the quadratures superposition Eq. (|24|) becomes classical. 
To understand the how the mode dynamics changes whcn 
going above threshold, we recali that below threshold the 
trajectory a\{k c ,t) occupies a circular region centered in 
zero, in the phase space given by its real and imaginary 
part. Above threshold the mean amplitude has a macro- 
scopic value increasing with the pump intensity and the 
distribution of fluctuations in phase and intensity quadra- 
tures is rather different, as shown in Fig. [|a. A trajectory 
for a mode during a larger time would describe a circle in 
phase space. 

Thc origin of these large fluctuations in the phase qua- 
drature is well known in the theory of single transverse- 
mode non-degenerate parametric oscillators (NDOPO) JÏgJ. 
Due to the diffusion of the difference of the signal and 
idler phases, the above threshold solution is not stable, 
and "cannot be analyzed correctly by the assumption of 
small fluctuations and methods of linearization" [|46| . The 
situation is similar to that of the làser above threshold, for 
which a correct analysis is performed using intensity and 
phase variables and not lincarizing in the diffusing (phase) 
variable. In Ref. |Ï7j an exact steady-state Wigner func- 
tion is calculated in the single transverse-mode NDOPO 
by adiabatically climinating thc pump. The phase diffu- 
sion in a mode is evident in the radial symmetry of this 
distribution. 

The single transverse-mode NDOPO is equivalent to 
a three mode model, describing the extended DOPO near 
threshold through the relevant fields ai(±fc c )and ao(0) 
p8[. The stationary signal is 



ai(x,t) = ai(k c )e lkcX + ai(-k c )e 

. 9 + +»_ 



2Mfc c )|e* 



cos 



(29) 



with ai(±fc c ) — |ai(±fc c )|e í6 ' ± . We recognize the effect of 
the phases 8± in the near field. The sum of these phases 
fixes the global phase of the signal, locked to the pump, 
while the arbitrary phase difference fixes the spatial posi- 
tion of the stripe pattern. In continuous systems, where all 
modes are taken into account, the diffusion of the phase 
difference can be interpreted as the action of the Gold- 
stone mode |Ï6), that is neutrally stable, giving a continu- 
ous translation of the pattern. This is particularly evident 
in Fig. at threshold (E = 1). 

In a linear treatment below threshold, the fluctuations 
of {9+ + (9_) have zero average and are damped (with 
Aq — 0). Increasing the pump we always observe small 
fluctuations, but the average changes. Fig. ||b shows that 
the mean value {6 + + (9_) increases from its zero value 
with the distance above threshold. This is a non-linear ef- 
fect due to the feedback of the signal on the pump. Above 
threshold the average of the non-linear term in the pump 
equation is not zero, so the pump is no longer real and this 
induecs a phase rotation in the signal. Therefore it should 
be expected that the strongest quadrature squeezing will 
be eventually found for a local oscillator phase that de- 
pends on this phase rotation. Here, however, we will rc- 
strict our attcntion to thc non-classical features associated 
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with intensity correlations. In fact this corresponds to a 
measure of quadrature squeezing as they result from the 
interference of a local oscillator given by the mean signal 
field with the squeezed fluctuations of the same mode. 



5.1 Intensity correlations in stripe patterns 

The stripe pattern formed above threshold is due to the in- 
terference of signal beams with opposite critical wavevec- 
tors. Momentum conservation leads to the entanglement 
between these signal beams [{7,14,21]]. This gives non-clas- 
sical intensity correlations characterized by V = —0.5. 

Incrcasing the pump intensity we observe excitation of 
harmonies of the critical wavenumber (compare far field 
(FF) in Fig. for E = 1 and for E = 1.1). In Fig. @ we 
show the real part of the near field (NF) pattern in the 
pump and in the signal and the corresponding FF inten- 
sities for pumps E = 1.02 and E = 1.1. We observe that 
the odd harmonies are excited in the signal and the even 
ones in the pump mode, with an exponential decay of en- 
ergy at higher wavevector modes. The presence of a multi- 
mode interaction means that the momentum conservation 
no longer constrains the intensities of the twin beams, 
as it does below and at threshold. However, as shown 
in Ref. plfl , we do observe in this regime the symmetry 
(-ZVi(fc)) = (Ni(—k)) in intensity averages and quantum 
correlations between the critical modes survive, as shown 
in Fig. H for 1 < E < 1.1. The secondary process of up- 
conversion of pairs of signal photons +k c (or equivalently 
— k c ) to form pump photons +2fc c ( — 2fc c ) does not seem 
to destroy the quantum correlations between the signal 
"twin" photons (+k c and — k c ). In principle this process 
gives an incoherent depletion of the "twin" beams, but 
probably due to the smallness of this secondary effect, 
the quantum correlations associated with the fundamen- 
tal process survives. 

The spatial spectrum of the intensity variance V is 
plotted in Fig. |8[ We observe that at 2% above threshold 
the spectrum is similar to the spectrum below threshold 
(compare with Fig. A reduction of the squeezing is 
observed (peak at k = 3fc c = 0.9 in Fig. ||), however, 
corresponding to the appearance of the third harmònic. 
Increasing the pump to 10% above threshold we observe 
an enhancement of the spectral bandwidth in which this 
reduction of squeezing appears. The third harmònic in the 
signal is involved in at least two important processes: the 
down-conversion of the homogeneous pump into twin pho- 
tons +3fc c and — 3k c and the secondary process of down- 
conversion of the second harmonies ±2fc c of the pump into 
opposite signal photons ±3fc c and ^k c . We note that as 
the Hamiltonian operator is Hcrmitian, the opposite (up- 
conversion) processes are also allowed. The observed re- 
duction of squeezing can be interpreted as a signature 
of the mentioned secondary process, in which the pairs 
signal +3fc c and — 3/c c photons are incoherently (not si- 
multaneously) generated and destroyed. In other words, 
signal modes are depleted independently, taking part in 
different cascading processes and gener ating harmonies. 
The entanglement should be preserved in opposite signal 



modes for which the fundamental down-conversion process 
prevails. We have confirmed that opposite spatial modes 
in the pump field do not show quantum correlations. 

We note that the variances V obtained with the Q- 
representation -after reordering- are in good agreement 
with the corresponding quantities calculated with the time 
dependent parametric approximation in the Wigner rep- 
resentation in Ref. Ml. 



5.2 Intensity correlations in spatially disordered 
structures 

With increasing pump intensity, a transition from a mod- 
ulated pattern to homogeneous solutions takes place |Ï9| . 
Due to the bistability in this regime, different homoge- 
neous solutions can be selected in separated spatial do- 
mains (see Fig. Q for E = 1.5). All solutions presented 
in Fig. [l] are obtained starting from a modulated initial 
condition at k c . Therefore we are stimulating , with a par- 
ticular initial configuration, the shape of the final struc- 
ture. In Fig. [9|a we show the stationary configuration for 
E = 1.5, in which the most excited mode is k c . The dis- 
ordered character of the structure gives rise to a broad 
spectrum, in which the other dominant modes are not har- 
monies of k c . The most intense signal modes combine to 
form pump modes: for example in Fig. [9|a the signal modes 
k c = AAk = 0.3 and 9Ak — 0.675 give the pump mode 
13Z\fc = 0.975 (50). The stationary disordered structure 
shows the same reflection symmetry (Ni(k)) = (Ni(—k)), 
of the stripe pattern considered in Sect. |5.l| . 

Studying the properties of the quantum fluctuations in 
this regime, we observe non-classical correlations between 
the two signal "twin" beams with critical wave-vectors, 
also in these spatially disordered structures (see Fig. || for 
1.2 < E < 1.5). This result does not depend on momen- 
tum conservation or on the presence of a regular pattern. 
We have also considered the entanglement prop erties of 
modes different from the critical one. In Fig. 10 we show 
the spatial spectrum of the variance V. As in the case of 
a regular stripe pattern, analyzed in the previous section, 
we continue to find quantum correlated twin beams. How- 
ever there are some interesting differences. The peak in V 
now corresponds to a strongly depleted (low mean inten- 
sity) signal mode for k = 0.975. This contrasts with the 
previous case, where k = 3k c was an excited mode. The 
most interesting feature is the appearance of a bandwidth 
of "twin" beams, where the signal field is intense; the cor- 
relations become classical for big wave-vectors (k > 1), 
where the signal is depleted more than the pump field, 
reaching asymptotically the level of coherent states (see 
dotted grey linc in Fig. |Ï0|). In conclusion, the "twin"- 
beams quantum correlations persist in disordered struc- 
tures. This signature of the fundamental down-conversion 
process is preserved throughout the region of intense sig- 
nal modes. 

The demanding question is how the spectrum of the 
variance V is influenced by the shape of the selected spa- 
tial structure. In this regime no special character is as- 
sociated with the critical wave-length periodicity. There- 
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fore we consider a stationary state of two domains, ob- 
tained from an initial step condition (Fig. |9)b) . Also in 
this case, with a very different stationary state, we ob- 
serve non-classical intensity correlations in the bandwidth 
< fe < 1, as shown in Fig. [ÏÏ^ (black line). In this case no 
peaks appear, suggesting that the presence of modes with 
reduced squeezing (as the peak for k = 0.975 observed 
in the previous structure) depends on the selected spatial 
structure. 

The key point is the relative importance of the funda- 
mental coherent process of twin photon down-conversion, 
and other incoherent cascading processes, which depends 
on the spatial configuration chosen by the system. To give 
morè evidence of this statement, we study the variance 

V in three different spatial configurations, obtained for 
E = 1.3, starting the simulations from noise, from rolls (as 
in Fig. 0a) and from a step function (as in Fig. ^b). In Fig. 
|ïl·| wc snow the spectral variances (a) and the mean FF 
intensity (b). The domain configurations obtained star- 
ing from noise and from a step function give overlapping 
smooth variances, while two peaks appear when starting 
from rolls . Also in this case the peaks appear associated 
with strongly depleted signal modes with low intensities, 
and can be in different positions depending on the selected 
spatial structure. The bandwidth of non-classical variance 

V seems to be a general feature, almost independent of 
the structure selected. It corresponds to the FF region of 
intense tilted signal beams. 



6 Conclusions 

Non-linear optical systems present a wealth of physical 
phcnomena including self-organizing spatial patterns and 
quantum correlations. Commonly employed methods to 
study these phenomena include the use of the positive P 
representation and the Wigner representation with third 
order derivatives neglected. Here we have investigated the 
use of the Q representation for studying quantum cor- 
relations in the DOPO at and above threshold. Positive 
diffusion is not guaranteed with the Q representation and 
this can lead to problems with divergent trajectòries. For 
the DOPO, however, only positive diffusion occurs unlcss 
the fluctuations are strong enough to push the pump field 
up to twice the threshold value. This never occurred in our 
simulations and we have not attempted to calculate the ef- 
fects of such highly unlikely trajectòries on our ensembles 

E- 

Below threshold entanglement in quadratures and non- 
classical intensity correlations are obtained, in agreement 
with linear analytical results. At threshold, a pair of quan- 
tum correlated twin beams is generated. These beams have 
wavevectors ±fc c , corresponding to the critical wavelength. 
The quantum correlations are a natural consequence of the 
fundamcntal microscopic process in which a single pump 
photon is converted into a pair of signal photons. As we 
move further above threshold the twin beams can recom- 
binc to generate new pump photons. The combination of 
a +k c and a — k c photon regenerates one of the original 



pump photons. The combination of a pair of +k c pho- 
tons or a pair of — k c photons, however, is a new pro- 
cess and introduces higher harmonies in the pump and 
thence in the the signal. Such processes can also degrade, 
but not completely suppress, the quantum correlations in 
some signal modes. This is a signature of the incoherent 
depletion of those modes. Yet further above threshold we 
enter a regime of spatially disordered structures. Remark- 
ably, quantum correlations persist even in this regime, in 
the bandwidth of intense signal modes, where they take 
a form that depends on the spatial pattern that is gener- 
ated. 
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Fig. 2. Normal ordered variance of the damped quadrature 
X-(k c ) normalized to shot noise: diamonds are results ob- 
tained by numerical simulation, while the continuous line cor- 
responds to the analytical expression Eq. (|2^). For any trajec- 
tory at given pump intensity, we average during a time of 10 7 , 
integratingwith a time discretization of 10 -3 (with time scaled 
as in Eq. (Q). 
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Fig. 3. Variance of the undamped quadrature X+(fc c ): the 
diamonds are results obtained with numerical simulation, while 
the continuous line corresponds to the analytical expression 
Eq. (]27|). At the last point, corresponding to E = 1, the linear 
treatment gives an infinite variance (the asymptotic behavior 
is represented by a dashed line) , while our non- linear treatment 
gives the expected saturation. 
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Fig. 4. Analytical (continuous line) and numerical (diamonds) 
twin beams correlations V(k) (Eq. (|Í|)) below threshold (E = 
0.99). The insert shows the mean intensity of the signal, with 
a maximum at k c = 0.3. 
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Fig. 5. Variance V(fc c ) at the critical wave- vector (k c — 0.3), 
increasing the pump intensity from below to above threshold 
(E c = 1). The diamonds are numerical results and the contin- 
uous line is the analytical value obtained by linearizing below 
threshold. 
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Fig. 6. a) Trajectory in the phase space of ai(k c ) during 10 7 
scaled units for E = 1.02. b) Phase sum < 0+ + 6- > (in 
radians) increasing the pump intensity. 



NF 



FF 



1. 1 




V V V V vvv w 



10 



2 







10 

-2 




10 



10 



10°'- 



-2 : 

-4 - 



í INI i 



41 
x 



82 



-1.5 



0.0 

k 



1.5 



Fig. 7. Snapshot of the real part of the near fïelds (NF) and 
intensity (log scale) of the far fields (FF) for both the pump 
and the signal. The grey plots are obtained for a pump value 
E — 1.02 and the black ones (in continuous or dotted lines) for 
E = 1.1. 
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Fig. 1. Spatiotemporal evolution of the real part of the near field (NF) and of the intensity of the far field (FF) (in log scale), for 
different vàlues of the pump E = 0.999, 1, 1.1, 1.5. The FF intensity is defined as \ai(k)\ 2 , where a;(fc) is the Fourier transform 
of the near field Qí(x). The horizontal coordinate is the transversal position (x in NF and k in FF) describes by 64 points, and 
the vertical one is the time interval 10 7 (in 7 units), using a discretization time step of At = 0.01. The initial condition for the 
signal is ai(a;,0) = 10 _5 (e(a;) + 10sin(fc c x)) with e{x) Gaussian random numbers of variance one. 
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Fig. 8. Twin beams correlations V(k) above threshold, for 
E = 1.02 (grey) and for E = 1.1 (black). 




Fig. 9. Snapshots of the real part of the NF and intensity (log 
scale) of the FF for the pump (thin line) and the signal (thick 
line), for E = 1.5, starting from a rolls pattern (a) and from a 
step function with vàlues — 1 and +l(b). 
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Fig. 10. Twin beams correlations V(k) for pump E = 1.5. The 
grey (black) line is the spectrum V(k) for the pattern shown 
in Fig. ^a (Fig. ^b ). The grey dotted line is obtained consid- 
ering a system with the same size (L = 4A C ) but with a fïner 
discretization (128 instead of 64 points), starting from a stripe 
configuration of critical wave-length. In this way we can see the 
asymptotic behavior of the spectrum for large wavevectors. For 
small wave-vectors the results are in good agreement with the 
simulation using 64 points. 
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Fig. 11. a) Twin beams correlations V(k) for pump E — 1.3, 
corresponding to different spatial structures, obtained starting 
from a step function (black line), from noise (grey line) and 
from a stripe pattern with critical wave-length (dotted line). 
Mean FF intensity in the pump (< No >) (b) and in the signal 
(< Ni >) flelds (c), for the same three spatial configurations. 



